Theory of the Bloch OsciUating Transistor 
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The Bloch osciUating transistor (BOT) is a device, where single electron current through a normal 
tunnel junction can be used to enhance Cooper pair current in a mesoscopic Josephson junction lead- 
ing to signal amplification. In this paper we develop a theory, where the BOT dynamics is described 
as a two-level system. The theory is used to predict current-voltage characteristics and small-signal 
response. Transition from stable operation into hysteretic regime is studied. By identifying the 
two-level switching noise as the main source of fluctuations, the expressions for equivalent noise 
sources and the noise temperature are derived. The validity of the model is tested by comparing 
the results with simulations. 

PACS numbers: 74.78.Na, 85.25.Am, 85.35.Gv 
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I. INTRODUCTION 

The Bloch oscillating transistor (BOT)i~— is a device 
based on tuning the probabiUty of interlevel switching 
in a mesoscopic Josephson Junction (JJ). The equiva- 
lent circuit is shown in Fig. 1(a). The current Ic at 
the collector(C) -emitter(E) -circuit is controlled by the 
base current Ib leading to transistor-like operation. The 
physics behind this is based on controlling the state of 
the JJ by means of quasiparticles tunneling through the 
normal tunnel junction connected to the base electrode 
(B). 
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FIG. 1: (a) Schematic circuit of a BOT connected to a source 
and a load. Source Rs is connected to the base electrode 
B and load Rl to the collector electrode C. The lead capaci- 
tances from the electrodes to the ground are Cb and Cc- The 
BOT itself consists of a Josephson Junction (JJ) between the 
emitter electrode E and the superconducting "island", normal 
tunnel junction between B and the island and large R ^ Rq 
resistor between C and the island, (b) The state diagram of 
the JJ and possible transitions induced by charge tunneling 
through the junctions. 



The state diagram as function of the "island" 
(quasi)charge Qi is shown in Fig. l(b)^, where also 
the transitions are illustrated. It is assumed that the 
Josephson coupling energy Ej is smaller or of the same 
order as the charging energy Ec — e^/2Cs, and that 
R, Rti , Rt2 ^ Rq ■ Here Cs is the total capacitance of 
the island, R is the collector resistor, and ii!Tiand Rt2 
are the tunnel resistances of the two junctions. The quan- 
tum resistance Rq = /i/4e^ « 6.5 kfi. When C is biased 
at a point, where Vc e/Cs, the charge tends to relax 
through the collector resistor Rc towards VcCs- Here 
Vc is the collector voltage. At Qi — e it is likely that 
a Cooper pair (CP) tunneling through the JJ returns 
the island back to Qj = — e. Repeating this cycle, the 
Bloch Oscillation^., leads to a net current through the 
CE-circuit. A competing process is the Zener tunneling-, 
which takes the system to the upper bands {\Qi\ > e 
in the extended band picture we are using). There the 
Cooper pair tunneling is blocked. One or more quasipar- 
ticles tunneling through the base junction may return 
the system back to the lowest energy band i.e. back to 
\Qi\<e. 

The Zener tunneling probability between 
bands n and n — 1 is given as -P^_i^„ ~ 

exp (- (7r/8) {Ej/Ec f (eEc/hlo)) . In the 

limit of small Ej it approaches unity very rapidly as 
n increases. Therefore we can neglect Cooper pair 
tunneling at higher bands, and consider the device as 
a two-level system. Below we will call the lowest band 
with n = the "first level" and higher bands with 
71 > 1 the "second level". Tunneling events, which cause 
transitions between the levels will be called interlevel 
transitions (or upwards and downwards transitions) and 
those, which cause transitions between higher bands, 
will be called intralevel transitions. Controlling the 
probability of downwards transitions leads to control of 
average Ic and thus to transistor- like characteristics. 

The BOT was recently experimentally realized^, and 
simulations showed that its properties can be quantita- 
tively predicted with a computational model^. It is po- 
tentially useful in cryogenic applications such as readout 
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circuits of radiation detectors, or measurement of small 
currents in quantum metrology. The aim of this article is 
to gain more insight into the dynamics of the BOX and to 
find the noise properties. To be able to do so, we derive 
an analytic theory, and study its applicability by com- 
paring the results to computational data. Small-signal 
parameters of BOT as well as equivalent noise sources 
and the noise temperature are derived. 



II. COMPUTATIONAL MODEL 

Since we assume that Ej <^ Ec we can interpret the 
quasicharge as the real charge, and the energy diagram 
will reduce to the simple parabola (dashed line in Fig. 
lb). We can account for the Josephson coupling per- 
turbatively and calculate the CP tunneling probabilities 
using the P (i?)-theorji2iS. This takes the electromag- 
netic environment into account. The consequence is that 
the tunneling does not happen strictly at Qj — e, but is 
rather represented with a finite distribution^. 

The simulation is done by integrating Qj in the time 
domain. The charge Qj — Q2 ~ Qi becomes 



Ib- In this case Vb is integrated as 



dVB 



dt 




(4) 



QPl, 



where Ib is the base current and Cb is the cable capaci- 
tance. Here we have assumed that Cb ^ Cs- 

Both Cooper pair and quasiparticle tunneling proba- 
bilities are obtained from the P (i?)-theory-, so that the 
effect of electromagnetic environment (the collector re- 
sistor) is also included. The given tunneling rates are 
towards the island. By reversing the signs of voltages 
Vi and V2 the rates from the island are obtained. The 
Cooper pair tunneling rate is 



Tcp {V2) 
and P (E) is defined as 



ttEI 



Pi2eV2)., 



(5) 



I dtexp (j it) + lEt], (6) 



dQi Vc - V2 



dt 



Rn 



dQi 
dt 



QPl 



dQi 
dt 



QP2 



dQi 
dt 



CP 

(1) 

where the first term is the collector current, the sec- 
ond and the third terms describe quasiparticle tunneling 
through the two junctions and the last term describes 
the Cooper pair tunneling through the Josephson junc- 
tion. The corresponding tunneling rates as the function 
of the state of the system are derived below. The volt- 
age V2 is across the JJ. The charges and voltages across 
both junctions, Qi (base junction) and Q2 (JJ) needed 
for the tunneling rate computation below are obtained as 
follows: 



Oi = 7^ {C1C2VB ~ CiQi) 

Q2 = 7^ {C1C2VB + C2Q1) 



(2) 
(3) 



The voltages across the junctions are Vi — Qi/Ci. In 
numerical simulations we will assume zero load resistance 
{Rl = 0), i.e. the collector voltage Vc = Vcc is fixed. 
In some simulations the base electrode is assumed volt- 
age biased (i.e. i?5 = and Vb = Vbb)- Note that 
below we will also use parameter Vg — Vb — Vc as base 
bias parameter for convenience. From the circuit point of 
view the choice between Vb and Vb is the detailed bias 
arrangement. 

In some simulations Rs will be assumed much higher 
than the input impedance. This fixes the base current 



where the phase correlation function J (t) is 



J (i) = 2 



duj Rc/Rqi 



uj 1 



^RcCy.) 
2kBTj 



(7) 



coth ( ^ .^"^ I (cos (ujt) — 1) — i sin (ujt) 



Here Rq2 = h/Ae^ sa 6.5 kil is the quantum resistance 
for Cooper pairs. For the quasiparticle tunneling formu- 
las given below value Rqi — h/e^ k, 25.9 kfi is used 
instead. Capacitance Ce ~ Ci + C2 is the total capaci- 
tance of both junctions, T is the temperature and fc^ the 
Bolzmann constant. 

The quasiparticle tunneling rates through the junc- 
tions [i = 1 for the base junction and i ~ 2 for the 
JJ) are 



Tqp. (V) = 



e^RT^ 



^ dEdE'^^^^i^^^ 
7Vi(0) N2{Q) 



(8) 



x/(i?) [l-J{E' + eV)]P{E-E'), 



where Rti are the tunneling resistances. The den- 
sities of states on both sides of the junctions are 
Ni (E) and N2 (E). For both NIS junction and the J J 
Ni (E) /Ni (0) = e (£;2 _ A2) E/VE'^ - A2. For the JJ 
N2 {E) IN2 (0) = e - A2) E|^/E^^ - A2 and for the 
NIS Junction N2 (E) /N2 (0) = 1. The superconducting 
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gap is A, / {E) is the Fermi function and 8 {E) is the 
step function. 

One should notice that the voltages appearing in the 
tunnel rate formulas are time dependent as opposed to 
the standard P (i?)-theory, where the fixed bias voltage 
is used. By doing so we can account for the fact that the 
tunneling probabilities depend on the state of the system. 



III. ANALYTIC THEORY 

In the theory derived below, BOT is modelled as a 
mapping of voltages Vb and Vc into currents Ic and Ib ■ 
We assume that a single tunneling event will not affect 
the voltages. This is the case, since Cb, Cc ^ in a. 
practical experimental setup. 

We assume that 1 < Ec/kT < Rc/Rq and Ej < 
Ec, which means that the Cooper pair tunneling rate 
(Eq. lO) reduces to a delta spike centered at \Qi\ — e^. 
This recovers our interpretation of the two-level system 
with IQ/I < e as the first level and \Qi\ > e as the sec- 
ond level. We also assume that C2 ^ Ci and neglect 
quasiparticle tunneling through the J J. Below unneces- 
sary subscripts for capacitances and charges are dropped, 
i.e. C = C2 = Cs, Rt = Rti and Q = Qj = Q2. We 
analyze only the regime, where Vc > e/C and < 0, 
since this is interesting for the amplifier operation. 

The collector current is written as 



i/rt 



i/rt + i/r 



(9) 



where the first term is the Cooper pair current through 
the Josephson junction and Ib is the single electron cur- 
rent through the base electrode. The transition rates 
between the two levels are and F|. The "saturation 
current" , i.e. current through the J J at the first level, is 
Is ~ '2efB, where is the Bloch oscillations frequency. 
The base current is 



i/rt + i/rx 



(10) 



where {N^} is the number of electrons needed to induce 
a downwards transition. Here we have neglected the pos- 
sibility of single-electron tunneling, when the system is 
at the first level. This is justified, since typically the 
voltage \Vi\ is below the gap voltage in that case. The 
Eqs. lO and (|10|l give general IV characteristics for the 
BOT. 

Between tunneling events dQ/dt = {Vc^Q/C)/R. 
By integrating from Q = — e to Q = e, i.e. over one 
Bloch period one gets 



Is 



2e 
RC 



In 



[Vc/Vq + I 
\Vc/Vq-1 



(12) 



where we have defined Vq — e/C. 

The upwards tunneling rate (the Zener tunneling) can 
now be written asifi 



Is 



2e (N) 



(13) 



where 



{N) = exp 



RC 



e Vc/Vq - 1 



- 1 



(14) 



is the average number of Cooper pairs in one sequence 
of Bloch oscillations. One sequence here means the time 
between tunneling down to the first level and tunneling 
back to the second level. The Zener avalanche current is 
4 = TreEj/8hEc. 

The downwards tunneling at low temperatures and for 
large R is exclusively due to single electron tunneling 
through the base junction. To calculate (TVe) and F| we 
first derive an approximation for Tgpi. The most gen- 
eral form is obtained from Eq. For many purposes, 
however, a piecewise linear approximation is sufficient: 



QPl 



1 



Tqpi (t) 



RtC 
1 

RtC 
0, if \Vi 



Qit) 



Vc 



Vq 



(15) 



if I'^^il > Vg 



where the gap- voltage is Vq, i.e., we neglect the leakage 
current at the subgap voltages. In the most straightfor- 
ward experimental realization the base junction is a NIS 
junction. In this case Vg = A + Vq/2 including the con- 
tribution of both superconducting and Coulomb gaps. In 
principle it is also possible to realize the base junction as a 
NIN junction, i.e. with suppressed superconductivity on 
the other electrode. Then the gap voltage is Vc — Vq/2. 
In general, Fgpi is a function of time due to the time 
dependency of charge. 

We now have to separate two different regimes to find 
analytic approximations for {N^) and F|. If Vc < 2Vq 
one electron always suffices to return the system to the 
first level, i.e., {Ne) = 1. In this case the probability 
distribution of the first quasiparticle tunneling event after 
the Zener tunneling is 



RC 



In 



Vc/Vq + 1 
Vc/Vq - 1 



(11) 



t/At 



P it) = Fqpi it) \im^ n (1 - Tqpi (jAt)) , (16) 
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which is the probabiHty that an electron will tunnel at 
time t times the probability that it has not tunneled at 
earlier times. The charge in Eq. H15|l before the first 
quasiparticle event obeys simple i?C-relaxation, i.e. 



Q (t) = CVc{l 



^-1 
Vc 



(17) 



where we assumed that Zener tunneling occurs at t = 
(or equivalently that Q (0) = e). The average rate is 
the inverse of the weight of the distribution given by Eq. 
(O. Thus 



(Ne) = 1 



r, = 



J,°°tP{t)df 



(18) 
(19) 



In general, the downwards rate has to be evaluated 
numerically from H19|) . However, if we further assume 
that the transient in Eq. ((T7|) is short, Eq. reduces 
to a simple exponential distribution. This is equivalent 
to assuming that Q (t) « CVc at all times, and thus 
Eq. I|15(l also becomes time independent. Now simply 
— Fgpi, and it follows 



(2.5i?T/i?+ 1-1) (V^/Vq)^ term in Eq. ^ dominate 
at Vc « 2Vq and large \Vq\. In this case only one 
quasiparticle is needed to induce a downwards transi- 
tion. This is possible, if the tunneling occurs during the 
transient immediately after the Zener tunneling, while 
still Q (t) < 2e. The terms are actually an approx- 
imation of Eqs. lO and lO at Vb « 2Vq. The 
exp(0.3 exp(...))-term dominates, when several tunneling 
events are needed to induce an interlevel transition. The 
very strong dependence is roughly explained as follows. 
Let us assume that 2Vq < Vc < 3Vq and the island 
charge is initially Q « CVc (see Fig. 1). Now at least 
two quasiparticles tunneling rapidly one after another are 
needed to induce a downwards transition. The quasipar- 
ticle tunneling probability according to Eq. IjlSfl is at 
its maximum, when Q w CVc- However, after the first 
tunneling event Q drops down to CVc ~ e and therefore 
the probability also drops. Hence the probability for the 
second quasiparticle to tunnel before the charge relaxes 
back to Q > 2e is small. The charge therefore tends 
to oscillate between Q k, CVc and Q « CVc — e for a 
long time before the rather improbable event at Q < 2e 
happens. This generates a large quantity of intralevel 
transitions thus increases (A^e) and F|. 



i 



CRt \Vq 2 



(20) 



If Vc > 2Vq it is possible that the first electron tunnel- 
ing through the base junction does not cause a transition 
to the first level, but some of the single quasiparticle 
events lead to intralevel transitions instead. This was 
found to have a dramatic effect in the experiments^ij^, 
and is found to be an important issue from the device 
optimization point of view as well. To solve {Ne) and F^ 
analytically from Eqs. (|15|l and H17|) in this case is unfor- 
tunately impossible. To find a sufficient approximation 
for our purposes, we have solved the problem numerically 
and searched for a proper fitting function. For simplicity 
we have assumed an NIN junction at the base electrode, 
i.e. Vc = Vq/2 in Eq. H15|l . Some fits are shown in Fig. 
2 and the result is 



0.04 



Rt 
~R 



(21) 



X exp 0.3 exp 1.8 



Vc 
Vo 



0.27 



VcVb 
V^ 




(22) 



The fit is accurate, when Rt < R- The weaker de- 
pendence indicated by the unity term in Eq. (|21|l and 



IV. COMPARING NUMERIC AND ANALYTIC 
IV CURVES 

In Fig. 3(a) we show a simulated set of Ic — Vc curves 
(open circles), where the base is voltage biased. The base 
voltage Vg is varied, while other parameters are i? = 1.5 
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FIG. 2: Computed data for (A^e) and obtained by solving 
Eqs. I15II and 11171 numerically (open circles). In each frame 
the base voltage is varied as Vg = — l.Oe/C, -1.5e/C, -2.0e/C, 
-2.5e/C, -3.0e/C from left to right. The lines correspond the 
fits, i.e. Eqs. ll?Tl and CT. 
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FIG. 3: (a) Computed Ic - Vc plots with R = 1.5 MQ., C = 
0.2 £F, Rt = 12 MSI, Ej/Ec = 0.1, T = 40 mK and A = 1.5 
mV (open circles) . The base voltage has been varied as Vg = 
-2.5e/C, -3.0e/C, -3.5e/C, -4.0e/C, -4.5e/C from down to 
top. Solid lines represent analytic values calculated from Eq. 
@ together with approximations from Eqs. 11811 and 1191 . 
Dashed lines are corrected analytic curves, which take base 
junction nonlinearity at the finite temperature into account, 
(b) Computed Ic — Ib plots for the same device (solid circles) 
at Vc =1.25e/C, 1.5e/C, 1.75e/C from up to down. The 
open squares shows the same simulation without " Cooper pair 
back-tunneling" and lines show analytic predictions. 

Mn, C = 0.2 fF, Rt = 12 Ej/Ec = 0.1, T = 40 
mK and A = 1.5 mV. Corresponding analytic curves 
(solid lines) are calculated from Eq. © using the ap- 
proximation of Eqs. IjlSI) and H19() when calculating (iVe) 
and Fj. The agreement is reasonably good. An error is 
introduced due to the assumptions made on the shape 
of Tqpi (Vi). If one uses the full temperature-dependent 
form (Eq. ©) instead of Eq. (|15() the agreement is im- 
proved especially at low values of as denoted by the 
dashed hne in Fig. 3. 



The simulated curves (solid circles) fall below the theo- 
retical curves (lines) Ic = (2 (N) + 1) Ib (see also Sec- 
tion VI), i.e. the current gain is suppressed. However, if 
we artificially forbid the " Cooper-pair back-tunneling" in 
the simulation (open squares in Fig. 3(b)) the agreement 
is clearly improved. This shows that the effect indeed is 
the main factor suppressing the current gain in the mode 
of operation governed by approximation given in Eqs. 
(|18ll and (|19|l . Another mechanism due to spontaneous 
downwards transitions was discussed in Ref^, but it was 
found to be insignificant in this case. 



As the tunnel resistance of the base electrode is de- 
creased and Josephson coupling increased in simulations 
and experimental^)^ it has been found that the active bias 
region moves towards higher Vc indicating that the ap- 
proximation of (iVe) and F^ given in Eq:s H21|) and 1)22(1 
becomes relevant. In Fig. 4(a) a set of simulations with 
parameters similar to those considered above, with ex- 
ceptions Rt = 375 kn, Ej/Ec = 0.2 and A = for 
the base junction (i.e. we have assumed that the base 
junction is a NIN junction here). At the upper set it 
is again shown a set of simulated and analytic Ic — Vc 
curves showing a reasonable agreement. The agreement 
is again further improved by forbidding the " Cooper-pair 
back-tunneling" in the simulation, which is shown in the 
lower set of curves. 

The devices analyzed above have relatively small C, 
which makes the voltages at which they are operated (of 




FJiaV] FJi.iV] 



The remaining disagreement can be found to be related 
to the temperature dependence of Cooper pair tunneling 
probabilities. Even if Ec/kT is as high as about 120, 
incoherent Cooper pair tunneling enhances Cooper pair 
current at Vc ~ Vq = 800 fiY. The lower value of simu- 
lated Ic at larger values of Vc was found to be due to the 
fact that after a Cooper pair tunnels to the island it can 
immediately tunnel out of the island due to incoherent 
Cooper pair tunneling. This effectively suppresses (N) , 
or equivalently enhances F|. The effect is especially visi- 
ble in Fig. 3(b), where a set of simulations with a current 
biased base electrode is performed for the same device. 



FIG. 4: (a) Computed Ic — Vc plots for a device otherwise 
similar to that of Fig. 3 except Rt = 375 kfl, Ej/Ec = 0.2 
and A = (an NIN base junction). The base voltages are 
Vg = -l.Oe/C, -1.5e/C, -2.0e/C, -2.5e/C from down to top 
(open circles) . Analytic IV curves (solid lines) are calcu- 
lated from together with approximations from Eqs. 12111 
and The upper set (lifted by 1.5 nA for clarity) shows 

the result with the full simulation model, while the lower set 
shows the result without "Cooper pair back-tunneling", (b) 
Analytic and computed Ic — Vc plots for a device having 
R = 500 kn, C = 1.2 fF, Rt = 250 k^, A = 400 ^V and 
Ej/Ec = 0.3. The two topmost sets have been lifted by 0.5 
nA and 0.8 nA for clarity. 
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order Vq) rather high. The advantage is that the temper- 
ature dependence is minimized as Ec/kT increases. The 
drawback is that higher band-gap materials are needed, 
since the voltage across the Josephson junction must be 
below 2A. The above devices could be realized using Nb 
technology, for which 2A 3 mV. More conventional Al- 
junctions have 2A « 400 /^V, whence capacitances have 
to be above or around 1 fF. In Fig. 4(b) a set of Ic — Vc 
curves are shown for a device with R — 500 kfi, C — 1.2 
fF, Rt = 250 Ml, Ej/Ec = 0.3. The topmost set con- 
sists of analytic curves, where at Vc < 2Vq w 270 /^V 
approximation of Eqs. (|18|l and (|19|l and at Vc ^ 2Vq 
approximation of Eqs. (|21|l and 122|l is used. The two 
lower sets are simulated at T = 20 mK and T = 300 mK. 
Although again qualitatively similar, at T = 20 mK the 
main source of disagreement is the enhancement of Fj at 
a finite temperature. At T = 300 mK the spike is spread, 
since at relatively large temperatures (now Ec / kT « 2.6) 
also Fj is increased due to incoherent Cooper pair tun- 
neling in a same sense as indicated in Refi^. 



T +U 



OVe 



(r 



BVb 



-e(iV) Jil^(^^- 

^^^^^^rt-fr/(iv,) dVc 

Fi 9Ft Ft 



aF| 



T^{T^+T^)dVc T^{T^+T^)dVc' 



Gqu 



FtF 



1 



dV 



Ft+Fx^Fi ir^+Ti)dVc 
1 9F| 1 dis . 



(28) 
(29) 
(30) 

(31) 



Tt (rt +Fx) dVc F| dVc 



c 



9x, 



where the fact that F|, {N) and Is are independent of 
Vb is used. Indices for constant quantities are dropped 
here for clarity. We have also defined 



V. SMALL SIGNAL PARAMETERS 

A complete small-signal model for the BOT can be 
j'iven as conductance matrix 



(23) 







Gout 


9m 




Vc 


is 




gx 


Gin 




VB _ 



Fx (Ft + Fx) f djNJ ^dVA 

In the approximation of Eq. 118|) Pb is zero, since (N^) is 
constant. Using Eqs. H21I) and (|22|l instead makes values 
/3b ^ ^ possible. We call (3b the "hysteresis parameter" 
of the BOT. 

For some purposes it is also useful to define current 
gain 



where icTiBiVCT^B are the small-signal components of 
collector and base currents and voltages, i.e. small vari- 
ations around the point of operation. The definitions of 
small-signal parameters are listed below: 



die 
dlf_ 



9m 

Gin 



LBJvc 

Using Eqs. (^5)) and this is given as 



(33) 



Gi: 



9x 

Gold 



diB 
dVB 

fdlc\ 
KoVb) 

diB 
dVc 



Vc 



Vc 



(dlc_\ 



(24) 
(25) 
(26) 
(27) 



Note that Vb is kept constant in the last two lines. This 
is the natural choice, if the circuit shown in Fig. 1(a) is 
used. However, if the emitter is voltage biased instead of 
the collector, V^ should be fixed instead. The choice does 
not have an effect on the analysis below, since we will 
mostly be assuming small whence Vc is constant. 
For completeness, the general formulas for small signal 
parameters are given here anyway. 

Using Eqs. lO and H10|) we can rewrite the parameters 



^ = i 



Is 



e Fx (N^) (1 - /3b) 



+ 1. 



(34) 



VI. AMPLIFIER PROPERTIES 

In this section we assume that Rl ^ i/\Gout\, i-e. 
that BOT is read out with a current amplifier and thus 
Vc is constant. For example, if one uses a dc SQUID 
as a postamplifier, Rl is close to zero. Even though a 
voltage amplifier is used i?^ w can be realized in prac- 
tice using current feedback. We will discuss two limits, 
one with approximation given in Eqs. 118|) and H20|) for 
evaluating Fx and (Ne). In this case Pb = 0. The second 
approximation uses Eqs. (|21|1 and H22(l to evaluate Fx 
and {Ne). Then it is possible to tune (3b close to unity. 
The emphasis is to find noise properties of the BOT at 
low frequencies. 

The noise current at the output of the BOT is obtained 
by assuming that the dominant noise mechanism is the 
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two-level switching noise. At low frequencies the corre- 
sponding sperctral noise density of the output current 
fluctuations is (see e.giii) 



The noise figure, defined as the ratio of total noise at 
the output divided by the noise contributed by the BOT, 



Si,out — 4/^ 



i-L T 



(35) 



where the collector current switches between values Ic = 
and Ic = Is- In other words we have neglected the 
single electron leakage from the base electrode here. 



The small-signal noise model for the BOT in the limit 
of small Rl is shown in Fig. 5. The signal and the 
noise from the source are described as current genera- 
tors isg and in,s in parallel with the source resistance 
Rs- The input and output impedances are i?i„ = 1/Gin 
and Rout = ^/Gout- The current generator (3iB at the 
output accounts for the gain. The noise added by the 
BOT is represented in a standard fashion (see e.gii^) by 
equivalent voltage and current noise generators (e„ and 
respectively) at the input. According to Fig. 5 the 
output noise of the BOT excluding the contribution of 
the source {in,s = 0) at the output is 



1/2 
i,out 



1 



R^n + Rs 



1/2 

en 



l/Rr. 



l/R^n + l/Ri 



1/2 



(36) 



where S'e„ and Sin are the spectral density functions cor- 
responding to e„ and z„, respectively. Note that e„ and 
in and are fully correlated with equal phases. Choosing 



5-1/2 



S 



1/2 




(37) 
(38) 



the generators become independent of Rs and produce 
the output noise in Eq. H35|l correctly. However, the 
backaction noise (i.e. the noise current iba through or 
voltage accross Rs) is not correctly predicted by the 
model. 




FIG. 5: A graphical representation of the small signal model 
of BOT in the limit of small Rl- The noise added by BOT is 
represented with equivalent sources in end Cn- 



f3Rs 



Rr. 



Rs 



Sin,. 



(39) 



where Sin,s — ^kTo/Rs is the spectal density function 
of in,s and To is a reference temperature. One gets opti- 
mum impedance Ropt and corresponding minimum noise 
temperature T„ by minimizing F and using the dcfiniton 
F = 1 + Tn/Ta. It follows 



R. 



opt 



Tn 



= \R,:. 



I Rin I 



(40) 
(41) 



The correlation of the two sources shows in Eq. H41|) in 
such a way that the prefactor is l/ks instead of l/2kB, 
which is the case for uncorrelated sources. The difference 
stems from the fact that now the amplitudes of the two 
sources rather than the powers are summed. 

If the approximation of Eqs. H18|l and (|20|l is used to 
evaluate {Ne) and Fj (whence also (3b =0), one gets for 
gain and noise parameters 



2(N) + 1 
2{N) + 1 



1/2 _ 



51/2 



R 



opt 



Tn 




(42) 
(43) 

(44) 

(45) 
(46) 
(47) 



In this mode the BOT acts as a simple "charge multi- 
plier" , where one electron trigs (N) Cooper pairs, thus 
(3 — 2 (N) + 1. The current noise can also be expressed as 
Sll^ = 2Ve7^(l + rx/r|)"\ InthelimitofsmaUFx/Fi 
the Bloch oscillation sequences are short compared the 
total length of the "duty cycle" l/Ti -f- l/r|. Then the 
equivalent current noise can be understood to be sim- 
ply the shot noise of the input current. In that case 



1/2 



The prefactor 2 instead of more famil- 



iar \/2 is due to the random length of charge pulses as 
opposed to standard shot noise. With large Fi/Fi, or 
with long Cooper pair sequences, the noise drops. The 
impedance also increases because single electron tunnel- 
ing is forbidden during the Bloch oscillations. One should 



8 



remember, however, that this is strictly true only in the 
absence of base junction leakage current. 

Lowering the equivalent current noise below the input 
current shot noise level can be understood as follows. 
In the limit governed by Eqs. (|42() - (|46|l the length of the 
duty cycle is determined by the base current Ib (Eq. I|1U|) 
with (A'e) = 1). If we have very short Cooper pair se- 
quences (or short l/T^) compared to the duty cycle, the 
output current is essentially a sequence of short charge 
pulses of size 2e (N) . This leads to shot type noise at 
the output, i.e. Slilt = 2^2e (N) Ic (Fig. 6(a)). The 
equivalent noise at the input is then Sj^^ ~ S',^{^j/2 (N). 

Using Ib w /c/2 (N) we then get S^jf^ = 2Ve7^- How- 
ever, if the base current is constant and the length of 
the Cooper pair sequence l/r^ (or equivalently {N)) is 
increased by increasing the Josephson coupling we even- 
tually have a situation, where l/Fj^ is the same as l/Tj 
was initially (see Fig. 6(b)). The current noise at the 
output is obviously the same is both cases, but the cur- 
rent gain 2 (TV) is larger in Fig. 6(b). Therefore the 
equivalent noise at the input must be smaller. 



As noted above, the spectral noise density of the back- 
action noise current {iba in Fig. 5) in general differs 
from Si,in- It can be shown, that for either T^/Ti <C 1 
or Fi/Fi ^ 1, it is exactly that of the base current shot 
noise, i.e. y/2elB- The maximum suppression of iBA 
occurs at F^ = F^, where the fano factor is 1/2. The 
reason for the difference in the equivalent current noise 
and the backaction noise is, that in the limit of large 
F|/F| the output current noise becomes fully anticorre- 
lated with iBA- Thus iBA does not directly determine 
the current resolution, or vice versa. To minimize the 
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FIG. 6: Collector current Ic schematically as function of time, 
when BOX is operated in the limit, where (Ne) = 1. In (a) 
Fi/Ft < 1 and (b) Ti/F^ > 1. The gray-shaded regions 
represent Cooper-pair sequences with area 2e (A'') . 



backaction noise, the device should be operated at a low 
base current. The low limit is here is set by sponta- 
neous downwards transitions due to incoherent Cooper 
pair tunnelingi. 

To investigate the parameters from the device opti- 
mization point of view, we assume a typical point of op- 
eration, where Vc = (3/2) Vq and = — Vq. Then it 
follows 




RT{l + {\n5)/2{RT/R) PY 



c.l/2 
5-1/2 



UE, 



. In5 R ^ 

, 1 H 

Rt \ 2 Rt 



'^EcRT{l+'f(^](3 



R. 



■opt 



Ec 
kB 



2 Rrj 



(48) 
(49) 

(50) 
(51) 
(52) 
(53) 



These formulas suggest, that Josephson coupling should 
be made large to maximize the current gain and to min- 
imize the added noise. One should also remember that 
the fluctuation effects mentioned in Section IV and Refi^ 
tend to suppress the gain which also increases the equiv- 
alent noise. Also the assumption Ej/Eq ^ 1 has to 
remain valid for the model to work. 

If the approximation from Eqs. H21|) and H22|l is used 
instead of Eqs. lfTH|l and for calculating F| and (A^e), 
the dominating terms are in many cases those dependent 
on 13b especially if /3b « 1. In Appendix A it is pre- 
sented a derivation of small signal and noise parameters. 
Approximations for F|, F^, {N^), (N), d (Ne) /OVb and 
dVi/dVB are made by eliminating the bias parameters at 
an interesting point of operation. "Fine tuning" of the 
device properties can be made by changing Pb, which in 
our approximation stands 



/3.^ 0.02 (A 



exp 




(54) 



while other quantities of interest are 
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FIG. 7: (a) and (b) show a computed set {R = 500 kf2, 
C = 1.2 £F, A = 400 AtV, T = 20 mK and RtIR = 0.5) of 
characteristic curves used to extract /3 and Ro-pt- The Joseph- 
son couphng is varied from Ej /Ec = 0.18 to Ej/Ec = 0.28 
(from right to left in (a) and down to up in (b). (c) Computed 
and analytic current gains (5 plotted against the optimum re- 
sistance Ropt as the Josephson coupling (or [3b) of the device 
is varied. The parameters are as above with the exceptions 
of varying Rt/R and T as shown in the legend. The bias 
point in the simulations with T = 20 mK is Vc = 3.5e/C and 
Vc = 4.5e/C for those with T = 300 mK. 



(55) 
(56) 

(57) 

(58) 

(59) 

(60) 

As (3b — 1 the current gain f3 diverges. However, the 
trade-ofF is that the optimum impedance Ropt also di- 
verges. The fluctuation at the output does not depend 

1/2 

on f3B, so the current noise and the noise tempera- 
ture T„ decrease at the same time. The transconductance 

1 /2 

gain gm and voltage noise Sen are independent of Pb- 




The physics in this limit can be understood as follows. 
With very large Pb the main effect of increasing Vb is in- 
creasing the number of electrons {Ne) needed to cause a 
downwards transition (see Eq. 1321) . This leads to decreas- 
ing Ib, i.e. negative input conductance. With very small 
Pb the only effect of increasing Vb is decreasing Tj . This 
leads to increasing Ib, i.e. positive input conductance^?!. 
At intermediate values, i.e. Pb ~ 1, the input conduc- 
tance is close to zero. The effect is that a small change 
in Ib causes a large change in Vb- Consequently F^, 
and thus also Ic change considerably. This leads to the 
enhancement of the current gain. 

A set of simulated Ic — Ib and Ib — Ve-plots with 
a varying Josephson coupling are shown in Figs. 7(a) 
and 7(b). The parameters were chosen so that the de- 
vice is realizable with Al-tunnel junctions (see the Cap- 
tion of Fig. 7). Current biased base electrode was as- 
sumed. This shows how the current gain and the in- 
put impedance increase without limit, as Pb approaches 
unity. As Pb exceeds unity the curves become hysteretic. 
If the source resistance Rs is large, hysteresis is a man- 
ifestation of negative input conductance. Therefore a 
sufficient stability criterion for all source resistances is 
Pb < ^- For small source resistances the device is sta- 
ble independently of Pb- The simulated IV curves be- 
come hysteretic at Ej /Ec ~ 0.25. According to Eq. 
Ej/Ec w 0.32 leads'to Pb = 1. 

We will next illustrate the trade-offs and test the va- 
lidity of Eqs. lfM|) - (|5n|l . The computed current gain as 
function of optimum resistance is shown together with 
the analytic approximation obtained from Eqs. H55|l and 
(|59|l in Fig. 7(c). Each of the three sets have different 
base resistance Rt and the bath temperature. Within 
each set Ej/Ec is varied. One can see, that the de- 
pendence P oc Ropt is correctly reproduced regardless of 
parameters, i.e. the property is quite generic. 

The current noise and the minimum noise temperature 
are shown as the function of the optimum resistance in 
Fig. 8. The computational noise data was obtained by 
performing a Fast Fourier Transform for the output cur- 
rent and averaging the low- frequency part, which gives 
Si, out- To further evaluate Sin, Ropt and T„, Eqs. H36|l . 
(|40|l and H41|l were applied with computed current gain 
P and input resistance Rin- Again the correct form of 
dependencies, i.e. Sl/^^ oc and T„ oc R~l are cor- 



rectly reproduced as compared to Eqs. and JHOJ). 

Differences in absolute levels can partially be explained 
through the inaccuracy of the numeric constants due to 
approximations made in Appendix A. To some extent the 
differences can also be understood with reference to ex- 
cess noise mechanisms discussed in Section VII. However, 
correct forms of dependencies and the order of magnitude 
are correctly predicted with the theory. 
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FIG. 8: (a) Computed and analytic results for the current 
noise spectral density Sin referred to input and (b) the min- 
imum noise temperature Tn as function of Ropt- The device 
and bias parameters are the same as in Fig. 7. 



VII. SUMMARY AND DISCUSSION 

We have developed an analytic theory to predict char- 
acteristic curves and noise properties of the Bloch Os- 
cillating Transistor. Even though it was derived at zero 
temperature, comparison to simulations at finite temper- 
atures showed reasonable agreement. The reason is that 
the main fluctuation mechanism, the two-level switching 
noise, is essentially temperature independent. Of other 
noise sources, for example, the thermal noise of the base 
junction is insignificant, since the junction is at a typical 
point of operation biased at eVi ^ kT. 

Two modes of operation were discussed. The small 
signal and noise parameters in both limits were obtained, 
i.e. Eqs. gHIl-lESl) and Eqs. (EB-lEOI)- In the first mode 
the device acts as a simple e— 2e (N) charge multiplier. In 
the second mode intralevel transitions play an important 
role. The consequence is the emergence of the hysteresis 
parameter f3B- This was found to have a drastic effect 
on device properties. It was shown that noise currents 
below 1 fA and noise temperatures below 100 mK can be 
obtained for optimum impedances of order a few Mfi. 

An additional noise mechanism is the shot-noise of 
the leakage current from the base-electrode. At a finite 
temperature also the bandwidth of the Bloch-oscillation 



increases^, and thus the shot-noise of the Cooper pairs 
adds to the total noise. This may explain the factor of 
about 3 increase in noise temperature as temperature is 
changed from 20 mK to 300 mK in Fig. 8(b). These 
effects appear at the output of the amplifier. Thus they 
add to the total output noise being additional terms to 
Eq. H35|l . Even in the presence of them the main conclu- 
sions of this article remain valid. 

We have also assumed arbitrarily large Cb, which is 
acceptable in low-frequency applications. However, if one 
wishes to increase the band, Cb should be decreased. 
The first effect of finite Cb is that the base voltage starts 
to fluctuate at frequencies typical to BOT dynamics. 

Most other well-known mesoscopic amplifiers, e.g. 
single-electron transistor (SET)ii or single Cooper pair 
transistor (SCPT)iS. are based on controlling current flow 
by charging a gate electrode. This is similar to the field 
effect transistor (FET), whereas BOT resembles a bipo- 
lar junction transistor (BJT) in the sense that a small 
base current is used to control a larger collector current. 
However, there are also important differences as well. For 
example, we have shown that the equivalent current noise 
of BOT can be brought below the shot noise of the in- 
put current. The reason is, that in BOT the noise at the 
output is partially correlated to the noise at the input. 

The l//-noise of the BOT is not addressed in this pa- 
per. However, as opposed to gate-controlled devices, the 
BOT is immune to background charge fluctuations. It 
is probable that the main 1 / /-noise mechanism is the 
fluctuation of the Josephson coupling. Due to symmetry 
considerations, it can be reduced by using bias reversal 
techniques typically used with a dc SQUID (— ). 
Authors wish to acknowledge fruitful discussions with P. 
Hakonen and R. Lindell. The work was supported by the 
Academy of Finland (project 103948). 



APPENDIX A: DERIVATION OF HYSTERESIS 
AND NOISE PARAMETERS 

In this Appendix we show the derivation of Eqs. H54|) - 
(|60|l . The aim is to derive an approximation, which ap- 
plies at /3b « 1 and Rt ^ R- 

The main task is to get sufficient estimates for Is, F^, 
F|, (Ne), d (Ne) /OVb and dTi/dVB near an interesting 
point of operation. We note that for the approximation 
from Eqs. (|2U and used for (A^'e) and F^ to apply 
the collector voltage must satisfy Vc ^ 2Vq. For simplic- 
ity we assume that Vc = 2Vq. When the system is at the 
second level, the voltage across the JJ is V2 ^ Vq and 
the collector current consists of the leakage current only, 
i.e. Ic = —Ib- Assuming that V2 = Vq and that the 
base junction roughly acts as a linear resistor we get by 
analyzing the circuit of Fig. 1(a) and noting that by def- 
inition Vb = V^+ Vc the result V^ = -Vq (1 + Rt/R), 
i.e. we have found estimates for the bias parameters. 
Using these and Eqs. H12|l . H13|) and H21|) we can readily 
write: 
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In the last Equation we have also utilized the assumption 
Rt < R- 

To find an estimate for we can use physical intuition 
and insight learned from simulations. Since the operation 
is based on switching between the two states, the system 
spends roughly as much time in both states. Thus 



Ti « r^. (A4) 

Although some error may be introduced by doing this, 
it is not too severe, since most properties depend more 
strongly on the derivative dTi/dVB, which will be calcu- 
lated separately. 

The derivative d {Ne) /OVb is obtained by direct differ- 
entiation of Eq. (|21|l . inserting the bias parameters Vb, 
Vc from above and applying the approximation Rt ^ R'. 
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To obtain dTi/dVB we first note that de- 
pends on Vb both through the explicit dependence 
in Eq. (|^ and through {Ng). Since {Ne) depends 
very strongly on Vb we will approximate dVi/dVB ~ 



{dTi/d{Ne)){d{Ne) /OVb). By differentiation of Eq. 
(|22|l . application of the bias parameters Vb and Vc from 

above and using Eq. (jsni 

we get 
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By using Eqs. (jA4|) and (|A2|) we now get 
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(A6) 



By inserting Eqs. IjAip - ljAGp into Eq. (|22Jl we get the 

estimate of the hysteresis parameter (3b , which is shown 
in Eq. l(5^ . 

The exponent (or (N)) in Eqs. HA2p and HA6|I affects 
on the device parameters mainly through Pb- For other 
purposes we may assume it roughly constant and solve 
it by setting /3b = 1 in Eq. (|^ . whence Eqs. lA2p and 
[6|) are simplified into 
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Eqs. lf5li|l - (|?^ are now calculated by substituting Eqs. 
(|X3|) . ((X4|) . ((X5|) . ((X7|) and l|X8)l into the definitions 
of interesting quantities, i.e. Eqs. (|23l, (jSHJ, 

(gOl and gH). 
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